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Abstract 

The structural, elastic and electronic properties of ReN are investigated by first-principles calcu- 
lations based on density functional theory. Two competing structures, i.e. CsCl-like and NiAs-like 
structures, are found and the most stable structure, NiAs-like, has a hexagonal symmetry which 
belongs to space group PQ^/mmc with a=2.7472 and c=5.8180 A. ReN with hexagonal symmetry 
is a metal ultra-incompressible solid and has less elastic anisotropy. The ultra-incompressibility of 
ReN is attributed to its high valence electron density and strong covalence bondings. Calculations 
of density of states and charge density distribution, together with Mulliken atomic population 
analysis, show that the bondings of ReN should be a mixture of metallic, covalent, and ionic bond- 
ings. Our results indicate that ReN can be used as a potential ultra-incompressible conductor. In 
particular, we obtain a superconducting transition temperature T c «4.8 K for ReN. 
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I. INTRODUCTION 



Transition metal nitrides are of intense interest for researchers.-^ 1 ^^ Not only do they 
offer the fundamental scientific challenge of finding methods to control, tune and enhance 
their peculiar physical properties, but they also have important technical applications such 
as cutting tools, 1,6 oxidation-resisting materials,- optical and magnetic apparatus.- Transi- 
tion metal nitrides possess particular mechanical, optical, electric, and magnetic properties, 
which attribute to their unusual electronic bonding (a mixture of covalent, metallic, and 
ionic bondings) and strong electron-phonon interaction. The higher valence electron density 
(VED) and the stronger hybridization between transition metal d electrons and nonmetal p 
electrons as well as less elastic anisotropy make some of them exhibit large incompressibility 
and high hardness, considered as powerful candidates of superhard materials. 8 

At present, there are three kinds of transition metal nitrides confirmed: X3N4 (X=Hf, 
Zr), 3 YN 2 (Y=Pt, Ir, Os),-^ and ZN (Z=Ti, Ta, Nb, Mo).- 6 -^ Many works, including deter- 
mination of the structure, discussion of elastic and electrical properties, have been carried 
out theoretically and experimentally.-^^^^^^^^'i 6 - The interaction of experiment and 
first-principles calculation has become a powerful tool to predict, test and confirm the struc- 
ture of the new materials.- For example, jointing experiment and first-principles calculation, 
Os^,- 1 ^^ 3 - 1 ^ IrN^^ and 5-MoN-^^ were identified as orthorhombic Pnnm, mono- 
clinic P2\C and hexagonal PQ^mc structures, respectively. The sites of nitrogen atoms in 
unit cell can be confirmed by theoretical calculation, which makes it possible to deeply study 
the mechanical, electronic and magnetic properties of transition metal nitrides. Haq and 
Meyer— reported the superconducting and electrical properties of Rei_ x N z (x=0.13~0.5) 
by means of ion implantation experimental technique. The rhenium nitride with fee struc- 
ture was observed and the corresponding superconducting temperature T c value of 4.5-5.0 K 
was obtained. However, Isaev et al 18 pointed out that transitional metal mononitrides with 
ten and more valence electrons in the NaCl-type structure are dynamically unstable which 
has been verified in MoN and NbN.— ^ As such, ReN with face centered cubic structure 
maybe unstable mechanically. So the crystal structure of ReN remains an open question to 
date. Moreover, to the best of our knowledge, the mechanical property of ReN has not been 
reported. 

In this paper, the structural, mechanical and electronic properties of ReN are investi- 
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gated by first-principles calculations. Because most of the transition-metal mono-nitrides 
are known to be generally based on cubic or hexagonal metal sublattice, and N atoms oc- 
cupy octahedral or trigonal prismatic interstitial sites, therefore, all the possible cubic and 
hexagonal structures are chosen as candidate structures of ReN, including zinc blende (ZB) 
(space group F43m), rocksalt (NaCl) (space group Fm3m), CsCl (space group Pm3m), WC 
(space group P6m2), and NiAs types (space group PQ^/mmc)M Then incompressibility and 
elastic anisotropy are discussed. Further, the superconduction transition temperature T c of 
ReN is estimated using the modified McMillan equation by Allen and Dynes.— Finally, the 
electronic property of ReN is given by electron structure calculation and Mulliken atomic 
population analysis. 

II. COMPUTATIONAL DETAILS 

Our calculations on structural and electronic properties of ReN are based on the full- 
potential linearized augmented plane waves plus local orbitals method (FP-LAPW) 2 - as 
implemented in the Wien2k package to solve the scalar-relativistic Kohn-Sham equations.— 
For the exchange correlation energy functional, the local density approximation (LDA) is 
employed. The muffin-tin radii are set to 1.9 and 1.60 bohr for Re and N, respectively. The 
maximum value l max for the wave function expansion inside the atomic spheres is limited to 
10. We expand the basis function up to RmtK max=~t '-0 (K max is the maximum modulus of 
the reciprocal lattice vector). Full relativistic approximation is used for the core electrons, 
and scalar relativistic approximation is used for the valence electrons of 5s5p4f5d6s for Re 
and of 2s2p for N. The self-consistent cycle is achieved by taking 1200 points in the first 
Brillouin zone. The convergence has been followed with respect to the energy and density. 

The elastic constants of ReN are obtained within the framework of the finite strain 
technique by CASTEP code using first-principles plane-wave basis pseudopotential method 
(PW-PP) based on DFT.— All the possible structures are also optimized by the BFGS 
algorithm 24 which provides a fast way of finding the lowest energy structure and supports 
cell optimization in the CASTEP code. In the calculation, the interaction between the 
ions and the electrons is described by using Vanderbilt's supersoft pseudopotential with the 
cutoff energy of 310 eV. In the geometrical optimization, all forces on atoms are converged 
to less than 0.002 eV/A, all the stress components are less than 0.02 GPa, and the tolerance 
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in self-consistent field (SCF) calculation is 5.0xl0~ 7 eV/atom. Relaxation of the internal 
degrees of freedom is allowed at each unit cell compression or expansion. From the full elastic 
constant tensor we determine the bulk modulus B and the shear modulus G according to 
the Voigt-Reuss-Hill (VRH) approximation.— The Young's modulus, E, and Possion's ratio, 
v, can be calculated by the formulae E = ^r G , v = ^b+g) • 



III. RESULTS AND DISCUSSIONS 



The equilibrium volume, bulk modulus and its pressure derivative are obtained by fitting 
the total energy calculated at different volumes to the 3-rd order Birch-Murnaghan equation 
of state.— The equilibrium structural parameters, density, valence electron density, the total 
energy, as well as the bulk modulus and its pressure derivative obtained from FP-LAPW 
method are listed in Table I. The energy- volume curve is plotted in Fig. 1. ReN is more 
energetically stable in the hexagonal structure under ambient pressure than cubic structure 
and NiAs-ReN structure is most stable, as shown in Table I and Fig. 1. In NiAs-ReN (space 
group P&z/mmc) the unit cell contains two chemical formula units (f.u.) in which two Re 
atoms are at 2a (0, 0, 0) and two N atoms at 2c (|, |, |) sites. The equilibrium volumes per 
formula unit (f.u.) are 23.5037 A 3 for ZB-ReN, 19.5258 A 3 for NaCl-ReN, 18.9065 A 3 for 
CsCl-ReN, 19.4775 A 3 for WC-ReN, and 19.0162 A 3 for NiAs-ReN, which are higher than 
that of Re metal (14.7173 A 3 )2£>2& due to the adding of N atoms. We note that the expansion 
of lattice when incorporating the N in the pure Re is not high enough to compensate for 
the decrease for the interatomic distance between the first neighbors. The nearest-neighbor 
distance d between atoms in ReN is much less than that of Re (which is about 2.74 A), 
which means that ReN maybe possess the much higher bulk modulus than pure Re metal 
(discussed later). The densities of ReN are 14.1476, 17.0299, 17.5877, 17.0721, and 17.4862 
g/cm 3 , for ZB-ReN, NaCl-ReN, CsCl-ReN, WC-ReN, and NiAs-ReN, respectively. 

The valence electron shell of Re is 5d 5 6s 2 , and that of N is 2s 2 p 3 so the total valence 
electron number is 12 per ReN molecule. The calculated VEDs are 0.5106 electrons/A 3 
for ZB-ReN, 0.6146 electrons/A 3 for NaCl-ReN, 0.6347 electrons/A 3 for CsCl-ReN, 0.6162 
electrons/A 3 for WC structure, and 0.6310 electrons/A 3 for NiAs-ReN. It is worthy to 
note that all VEDs for these three structures are higher than that of Re metal (0.4761 
electrons/ A 3 )^ and are compared with 0.70 electrons/A 3 for diamond.— Generally, the 
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higher VED, the stronger the incompressibility. In addition, from Table I, it is clear that 
the bulk modulus of ReN is larger than that (360 GPa) 2 - of Re metal except for ZB-ReN 
(unstable mechanically). Thus, we conclude that the ReN is less compressible than the 
pure Re. It is delighted that the estimated bulk modulus of ReN with NiAs, CsCl, and 
WC structures is even comparable to that of diamond, 30 which means that the ReN with 
these structures has a stronger incompressibility, indicating that ReN can be a candidate of 
superhard material. 

In order to get more accurate elastic constants of ReN using PW-PP scheme, all consid- 
ered structures are optimized firstly." The calculated elastic constants are presented in Table 
II. For a stable hexagonal structure, the five independent elastic constants c^- (en, c 33 ,cu,cn 
and C13) should satisfy the well known Born-Huang criteria for stability: 32 

C12 > 0, c 33 > 0, c u > C12, c 44 > 0, (c u + ci 2 )c 33 > 2c 2 3 ; 

while for cubic crystal, the three independent elastic constants (cn,ci2 and C44) satisfy 
inequalities: 32 

c 44 > 0, cu > |c 12 | , cu + 2ci2 > 0. 

Clearly, our calculated elastic constants Cjj for both cubic CsCl-ReN and the hexagonal 
NiAs-ReN satisfy the Born-Huang stability criteria, suggesting that they are mechanically 
stable. Cubic ZB-ReN and NaCl types are unstable mechanically because of their negative 
C44. Our results stand by for Isaev's viewpoint of which the transition metal mononitrides 
with ten and more valence electrons in the NaCl-type structure are unstable mechanically. 
WC-ReN type is also unstable mechanically due to the very small C44 value (near zero). In 
the following, only CsCl and NiAs types of ReN are discussed. 

From the elastic constants calculated above, bulk modulus B, shear modulus G, Young's 
modulus E, and Possion's ratio v are obtained and listed in Table II, which are important 
in order to understand the elastic properties of ReN. The bulk modulus of ReN calculated 
by using elastic constants agrees well with the one obtained through the fit to the 3rd order 
Birch-Murnaghan EOS, which shows that the discussion-above on compressibility of ReN 
is reasonable. To further discuss the incompressibility of ReN, the evolution of the volume 
compressions as a function of pressure are plotted in Fig. 2. For comparison, the volume 
compression of diamond is given simultaneously. It is evident that ReN with NiAs structure 
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is slight higher incompressible than diamond, while ReN with CsCl structure has slightly 
less incompressible than diamond in the entire pressure range (Fig. 2(b)). 

To explain the directional dependence of compression, fraction axis compression as a 
function of pressure for diamond and two structures of ReN are given in Fig. 3. Compression 
of the axes presents interesting anisotropy. For hexagonal NiAl-ReN, the c axis is more 
incompressible than the a axis. The incompressibility along a-axis is slightly less than that 
of diamond, while c axis is even more incompressible than the analogous axis of diamond. 
While for cubic CsCl-ReN, it is slightly lower incompressible than diamond (Fig. 2 and Fig. 
3). Therefore, we can claim that ReN is an ultra-incompressible solid. The shear modulli 
obtained are 248 GPa for CsCl- ReN and 237 GPa for NiAs-ReN, and corresponding G/B 
values are 0.588 and 0.526. The Possion's ratio is 0.25 and 0.28 for CsCl and NiAs types, 
respectively, which indicates that ReN has central interatomic forces and is relatively stable 
against shear.— The high G/B ratio or, equivalently, the low Poisson's ratio also implies a 
high value degree of covalency. 

It is well known that the elastic anisotropy of crystals has significant application in 
engineering science since it is highly correlated with both the hardness of materials and the 
possibility to induce microcracks in the materials. Hence it is important to calculate elastic 
anisotropy in order to understand the properties and improve the durability of materials. 
The anisotropy factor for cubic CsCl-ReN, v4=(2c 44 + Ci 2 )/cn=0.45, shows that cubic CsCl- 
ReN is elastic anisotropy. The anisotropy of linear compression along the a axis and c axis 
with respect to the b axis can be evaluated by means of two parameters: As a and ^4s c ,— 
where a value of 1 indicates compression isotropy and any departure from one corresponds to 
a degree of compression anisotropy. For low symmetry crystal, the percentage of anisotropy 
in compressibility and shear can be derived from two expressions, bulk modulus anisotropic 
factor Ab and shear anisotropy factor Ac, respectively— The calculated linear compression 
parameters A^ a , Ab c and two other anisotropy factors Ab and Aq are listed in Table III. 
NiAs-ReN is almost isotropy in compressibility (0.317%) and has slightly shear anisotropy 
(^4^=4.199%). From Table III, one can see that cubic CsCl-ReN is compression isotropy, 
and hexagonal ReN's c axis has much stronger incompressibility than that of a axis, which 
is in good agreement with fractional axis compression analysis mentioned above. The small 
elastic anisotropy of NiAs-ReN, including linearized compressibility and shear anisotropy, 
along with the ultra-high bulk modulus and the higher shear modulus (or c 44 ), indicates 
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that NiAs-ReN can be as a good candidate of superhard material. Comparing C44 with 
shear modulus for hexagonal NiAs-ReN, we find that they are almost the same, imaging 
that NiAs-ReN hardness can be given a consistent estimation by the viewpoints on hardness 
from Teter— and Jhi et alM. Here, using the correlation between shear modulus and Vickers 
harness given by Teter, the estimated Vickers hardness of ReN are 33.6 GPa for CsCl-ReN 
and 32.1 GPa for NiAs-ReN, and 36.2 GPa for MoN-ReN. All the hardness of ReN presented 
here are comparable to that of B^O (35±5GPa) and that of TiB2 (33±2GPa).— 

The Debye temperature correlates with many physical properties of solids, such as elastic 
constants, specific heat, melting temperature and superconduction transition temperature. 
Since the vibrational excitations mainly arise from acoustic vibrations at low temperature, 
the Debye temperature obtained from elastic constants is the same as that determined from 
specific heat measurements. The Debye temperature can be derived by bulk modulus, shear 
modulus and density p. 36 The calculated Debye temperature O^s are as high as 968 K 
for CsCl-ReN and 950 K for NiAs-ReN, which indicates that ReN maybe a superconduc- 
tor. So, we further estimate the superconducting transition temperature of ReN with NiAs 
structure using the modified McMillan equation,— where the electron-phonon coupling con- 
stant A is calculated by density-functional perturbation theory by means of the QUANTUM 
ESPRESSO package.— The calculated A is 1.67 and the estimated superconducting transi- 
tion temperature T c is about 4.8 K for the screened Coulomb pseudopotential //* equal to 
0.1, which agrees well with available experimental result.— 

To understand the correlation between the electronic properties and the mechanical prop- 
erties, we present density of states (DOS) and the electron density distribution of ReN in 
the equilibrium geometry. Electronic structure properties play an important role in material 
physical properties. The total density of states (DOS) and partial DOS of two phases of ReN 
are shown in Fig. 4 and Fig. 5. The obtained N(Ep) and the linear specific heat coefficient 
7 are given in Table I. The finite value N(Ep) shows that ReN is a metal. For two phases, 
the electrons from Re-d and the N-p states both contribute to the DOS near the Fermi level. 
The peak emerging in the lower energy region of the DOS curve mainly originates from the 
localized s states of N. The energy region just above Fermi level is dominated by unoccupied 
Re d states. The DOS of Re-d and N-p are energetically degenerate from the bottom of the 
valence band to the top of conduction band, indicating the possibility of covalent bonding 
between Re and N atoms. The covalent characteristic between Re and N atoms can be 
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confirmed by the charge density distribution. The charge density distribution in (110) plane 
for CsCl-ReN and (1120) plane for NiAs-ReN are shown in Fig. 6 (a) and (b), respectively. 
It is clearly seen that a strong directional bonding exists between Re and N in the ReN, 
which contributes to high bulk modulus B, namely, strong incompressibility of ReN. Ad- 
ditionally, the Mulliken atomic population analysis is carried out by using CASTEP code. 
The obtained total charge transfer from Re to N in CsCl-ReN (NiAs-ReN) is 0.55 (0.62), 
implying that the chemical bondings between Re and N have some characteristics of ionicity. 
Thus our results demonstrate that the bondings should be a mixture of covalent, mental, 
and ionic attribution in ReN. 



IV. CONCLUSION 



The structural, elastic, and electronic properties of different structures for ReN are in- 
vestigated based on first-principles calculation under the framework of density functional 
theory within local density approximation. Two competing structures, i.e., CsCl-ReN and 
NiAs-ReN, are found. High valence electron density, strong incompressibility and low elastic 
anisotropy indicate that ReN is an ultra-incompressible superhard solid. Electron structure 
calculation shows that ReN presents an obvious metal feature. In particular, we obtain a 
superconducting transition temperature T c ^4.8 K for hexagonal ReN. The strong hybridiza- 
tion between metal d-electron and nonmetal ^-electron are observed, indicating that there 
is a strong covalent bonding between Re and N. From Mulliken atomic population analysis 
along with discussion of mechanical and electronic properties, it can be concluded that the 
bondings in ReN should be a mixing of metal, ionic and covalent characteristics. Our results 
indicate that ReN can be used as a potential ultra-incompressible conductor. We hope this 
work can stimulate the experimental research for ReN. 
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FIGURE CAPTIONS 



Fig.l: (Color online) The total energy per ReN molecule as a function of volume (EOS). 

Fig. 2: (Color online) Volume as a function of pressure with respect to equilibrium volume. 
Vq is the equilibrium volume. 

Fig. 3: (Color online) Fractional axis compression as a functional of pressure. a , b , and c 
are lattice parameters of equilibrium volume at zero pressure. 

Fig. 4: (Color online) Partial and total DOS for the CsCl phase of ReN obtained from 
FP-LAPW calculations. Vertical dotted lines indicate the Fermi level. 

Fig. 5: (Color online) Partial and total DOS for the NiAs phase of ReN obtained from 
FP-LAPW calculations. Vertical dotted lines indicate the Fermi level. 

Fig. 6: The valence electron density contour for ReN obtained from FP-LAPW. Charge 
density is in an increment of 0.01 e/a.u. 3 from 0.01 e/a.u. 3 to 0.12 e/a.u. 3 . (a) (110) plane 
for CsCl-ReN; (b) (1120) plane for NiAs-ReN. 
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TABLE CAPTIONS 



TABLE I: Equilibrium lattice parameters, V (A 3 ), a (A), c (A), c/a, the shortest Re-N 
distance d, density p (g/cm 3 ), valence electron densities p e (electrons/A 3 ), bulk modulus B 
and its pressure derivative B' Q , relative total energy E tot (eV), DOS at Fermi level N(E F ) 
(States/eV/cell), and the linear specific heat coefficient 7 (mJ/mol.cell.K 2 ). V and E tot 
are of per chemical f.u. 

TABLE II: Zero-pressure elastic constants (GPa), the isotropic bulk modulus B (GPa), 
shear modulus G (GPa), Young's modulus E (GPa), Possion's ratio v, average elastic wave 
velocity v m (m/s), and Debye temperature 9^. 

TABLE III: The bulk modulus along the crystallographic axes a, b, and c (B a , B b , and 
B c ) for ReN. Percent elastic anisotropy for shear and bulk moduli A G (in %), A B (in %) and 
compressibility anisotropy factors A Ba and A Bc for ReN obtained from LDA calculations. 
Here, A Ba =f , A Bc =f. 
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TABLE I: 



Li et al. 
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TABLE II: Li et al. 



en c 33 c 44 C12 C13 B G £ v Q D 
ZB-ReN 297 -715 356 336 -255 -1025 - 
NaCl-ReN 456 -334 372 400 -27 -83 - 
CsCl-ReN 1015 164 124 421 248 622 0.25 968 
WC-ReN 640 654 0.55 240 368 429 53 152 - - 
NiAs-ReN 712 897 236 314 278 450 237 604 0.28 950 



TABLE III: Li et al. 



B a B b 


Be A Ba 


A Bc A B A g 


CsCl-ReN 1264 


1 


11.55 


NiAs-ReN 1237 


1627 1 


1.3155 0.317 4.199 
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FIG. 1: Li et al. 
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FIG. 2: Li et al. 
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FIG. 3: Li et al. 
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FIG. 4: Li et al. 
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FIG. 5: Li et al. 
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FIG. 6: Li et al. 
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